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We extend a treatment of the causal structure of space-time to active- 
sterile neutrino transformation-based schemes for lepton number generation 
in the early universe. We find that these causality considerations necessarily 
lead to the creation of spatial domains of lepton number with opposite signs. 
Lepton number gradients at the domain boundaries can open a new channel 
for MSW resonant production of sterile neutrinos. The enhanced sterile neu- 
trino production via this new channel allows considerable tightening of Big 
Bang Nucleosynthesis constraints on active-sterile neutrino mixing, including 
the proposed i/^ solution for the Super Kamiokande atmospheric 

deficit, and the four- neutrino schemes proposed to simultaneously fit current 
neutrino experimental results. 
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It is well known that resonant MSW (Mikheyev-Smirnov-Wolfenstein) |l|] transitions 
between active neutrinos and sterile neutrinos in the early universe could generate large 
lepton number asymmetries in the neutrino sector [@-|^. Here the lepton number for an 
active neutrino species Va is defined to be L^^ = {riu^ — np^)/n^, the net asymmetry in Ua 
over number density normalized by the photon number density n^. 

Not well noted is a crucial feature involved in the lepton number generation process : 
that the lepton number asymmetry is first damped to essentially zero by the active-sterile 
neutrino mixing, then oscillates chaotically with an progressively larger amplitude as the 
mixing goes through resonances, until the asymmetry converges to a growing asymptotic 
value that is either positive or negative. As a result of this feature, the sign of the lepton 
number asymmetry is independent of the initial conditions which obtain before the instability 
begins, and is exponentially sensitive to the parameters involved during the chaotic oscillitory 
phase. In turn, the lepton number generated in this process may not have a uniform sign 
in different causal domains. Obviously, an upper bound to the size of these domains is 
the particle horizon H^^ ^ (90/87r'^)^/^5f^^/^mpi/T^ where g is the statistical weight in 
relativistic particles, T is the temperature of the universe, and rripi ^ 1.22 x 10^^ eV. 

The typical size of these domains is at least as large as the diffusion length of neutrinos 
at the time of the lepton number generation. These leptonic domains can persist as long 
as the resonant neutrino transition is capable of efficient lepton number generation. This 
is because any reduction of lepton number at domain boundaries due to mixing will be 
quickly reversed by the generation process so that a boundary region is incorporated into 
one domain or the other. 

The existence of leptonic domains in the early universe provide a new channel of produc- 
ing sterile neutrinos, via the resonant MSW conversion of active neutrinos to sterile neutrinos 
at domain boundaries, where lepton number gradients exist. If during the Big Bang Nu- 
cleosynthesis (BBN) epoch the population of sterile neutrinos from this channel becomes 
comparable to that of an equilibrated active neutrino fiavor, the extra neutrino degrees-of- 
freedom (an increase in g) can increase the primordial helium abundance significantly [§. 
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Therefore, such a production mechanism can be constrained by the observationally-inferred 
primordial hehum abundance. 

The details of the active-sterile neutrino transformation process in the early universe 
and the associated generation of lepton number asymmetries can be found in a number of 
previous works Here we only briefly summarize. In all calculations, we employ the 

natural units h = c = ks = I- A two-family system Ua ^ Ug {a = e, fi or r, and Us is the 
sterile neutrino) has a 2x2 evolution Hamiltonian H with Haa = Vz, Has = = Vx + iVy, 
Tiss = 0. The effective vector potential V during the BBN epoch (below the QCD phase 
transition temperature T < 100 MeV) is 

K = -^sin2^, V; = 0, K = -^cos2^ + ri + rJ, (1) 

where 5m? ml — ml , 6 is the vacuum mixing angle, and E is neutrino energy. The 



matter- antimatter asymmetry contribution to the effective potential is [10 



^ ±0.35GfT' [Lo + 2L,^ + L,^] (2) 

with the sign for z/^ and the "— " sign for z/^. The quantity Lq represents the contribution 
from the baryonic asymmetry and electron-positron asymmetry, i.e., Lq ~ 10^^°. The 
quantity L^,^ is the asymmetry in other active neutrino species up. For simplicity, and 
with no loss of generality, we will assume L,^^ = unless explicitly stated otherwise. The 
contribution to V from a thermal neutrino background is ||10| 



^ , -SOGlET^ for a = e; 
VI ^{ (3) 
-22GlET^ for a = fi, t. 

is the same for both z/^ and z/^. 

Resonances occur when the diagonal elements of the Hamiltonian Ti. are equal, i.e., 
Vz = 0. This is when z/^ and z/^ become degenerate in effective mass and maximally mixed 
in medium. Very simplistically, when Vz evolves through zero, z/^ and Us swap their flavors if 
the resonance is adiabatic (jdl^^/dtl <^ V^), and remain essentially unaltered if the resonance 
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is non-adiabatic (|dl^/dt| ^ V^). Since the resonance condition is energy dependent, only 
neutrinos with Ej-^s are are resonant at any given temperature. 

When = (or neghgibly small compared to \6'm'^\/2E), resonances can occur for 
the Ua ^ Vs system only if < (i.e., niy^ > m^^). For 6m'^ > 0, and at high enough 
temperatures for < 0, it has been shown that active-sterile neutrino transformation can 
damp Lq + 2L,y^ to zero very efficiently. The amplification of Lq + 2L^^ starts for 5m? < 
only as the temperature falls below the critical temperature 



This is also the temperature where a significant number of z/^ go through the resonance. 
The instability of Lq + 2Ly^ at Tc results from the non-linear characteristics of the MSW 
resonances, and the feedback effect they have on the lepton number asymmetry [Q. 

Due to the same non-linearity, Lq + 2Ly^ (or L^^ once Ly^ ^ Lq) oscillates chaotically 
near T^, with an exponentially increasing amplitude. This chaotic feature was first pointed 
out by Shi in numerical studies of monochromic neutrinos undergoing MSW resonance 0, 
and was later found to apply to systems with a distribution of energies as well 0. Eventually, 
at a temperature slightly below Tc, the oscillatory behavor abates, and L^^ settles into one 
of two fixed points (two T^^ power-laws) 



The sign of the emergent L^^, however, is independent of the initial Lq and L^^ above T^, and 
is exponentially sensitive to the evolution of the potential V during the oscillatory phase. 
The sign is therefore chaotic and uncorrelated across causal domains. 

In the power-law regime, the growth of Lj^^ is driven by the resonant conversion of to 
z/g (if Ly^ < 0) or Ua to (if L^,^ > 0) in the sector of the neutrino energy spectrum with 
Eres/T ^ O.Oel^mVleV^I \LyJ-\T/l MeV)-^ which is in general < 1 |],|. In this regime, 

quickly becomes negligible because it scales as T^. The power-law solution, Eq. (^, for 
\Lu^ \ is very stable, in that a \Lu^ \ significantly deviating from this solution will be quickly 
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Tc ^ 22 W/leV^ MeV. 



(4) 
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damped or amplified until it converges to Eq. (||). (A larger \LuJ implies a small Ej-cs/T, 
and therefore a less efficient generation of \Li^^\] a smaller \Lu^ \ implies a larger Et-^s/T and 
a more efficient generation of The rate of this convergence is din |L,y^|/dlnT > 4 in 

the power-law regime. 

As \Li^^\ increases, E-^^^s/T will slowly increase. This is because keeping up with the 
approximate growth requires progressively more Va or to be resonantly converted 
into sterile neutrinos. Eventually, as E^-cs/T sweeps through most of the z/^ energy spectrum, 
the growth of \L^^\ tapers off near its physical limit. This limit is \Lu^\ = 3/8, when the 
entire z/q, or z/^ population has been converted into sterile neutrinos. In turn, the lepton 
number generation process ceases, and any inhomogeneity of lepton number arising from 
the process begins to be smoothed out by neutrino diffusion. 

Because of the chaotic behavior of the sign of L,^^ , domains of lepton number with oppo- 
site signs are expected to form at the epoch when T ^ Tc. At the domain boundaries, mixing 
between different domains tends to reduce the asymmetries and increases the thickness of 
the boundary regions. However, the resonant neutrino mixing tends to maintain the solu- 
tion in Eq. and so narrows the boundary region. The thickness of the boundaries, very 
crudely, is thus the diffusion length of active neutrinos within the time in which \Ly^ \ grows 
by an e-folding, /A. Taking the z/q, collision rate Ty^ ~ G'^T^, we obtain a boundary 
thickness relative to the horizon scale ct 



for T > 1 MeV when active neutrinos are still diffusive. At temperatures T < 1 MeV, 
neutrinos decouple from the plasma and free-stream at the speed of light. In this case 
6d ~ ct. 

The existance of lepton domains results in a gradient of L,^^, and therefore a gradient 
in Vz- More importantly, the varying Vz at the domain boundaries satisfies the resonant 
condition for most z/^ crossing the boundaries. This is because [as Eq. (j^) shows] the power- 



law solutions within each domain result from resonant transitions of z/^ or z/^ with Ej-es -C T 
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]. As a result, at domain boundaries Ej-cs becomes larger. This is because i^res oc 
Most of and z/„ therefore undergo resonant transitions to sterile neutrinos at domain 
boundaries. 

This new channel of sterile neutrino production can populate a significant sterile neu- 
trino sea with a number density comparable to that of an equilibrated active neutrino flavor 
if (1) the resonant conversion at the boundary region is adiabatic; and (2) this new produc- 
tion channel for sterile neutrinos does not provide negative feedback to the lepton number 
asymmetry generation process, and so does not compromise the domain structure of lepton 
number. Here we demonstrate that both requirements can be satisfied. 

The adiabaticity condition at the resonances is 







dt 





dr 



HT 



dT 



(7) 



evaluated at = 0. The second term on the r.h.s. is of order HV^ . This term is always 
small compared to the first term as long as the leptonic domains exist. Since the spatial 
gradient is expected to be smooth across the boundary (whose thickness is determined by 
the diffusion/streaming process), we can employ the average \dVz/dr\ across the domain 
boundaries, ~ Q.1GfT^{L\^J - Li-J)/6d [where L^^ ~ -L^,;) satisfying Eq. (§)]. In turn, 
the adiabaticity condition Eq. becomes 
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(8) 



The production of sterile neutrinos via adiabatic conversion of at domain boundaries 
will not have a negative impact on the L^^ generation process and the domain structure of 
lepton number. Consider a domain boundary with a lepton asymmetry L[^^ on one side, 
and L[,~^ on the other. The L^^^ side has fewer z/q, than u^, and some from resonant 
i^a ~* conversions. The L^~J side has the opposite, with more z/q, than and some 
from resonant Va conversions. When neutrinos cross the boundary from the L^J^ side 

to the L^~J side, the resultant production of sterile neutrinos due to Ua i^s and Ua i^s 
has no bearing on the asymmetry of z/q on the side. The existence of sterile neutrinos 
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do not hinder the L^,^ generation process until the sterile neutrino population is comparable 
in numbers to the population [Q]. The resonant conversion of Us —>■ i^a due to the crossing, 
on the other hand, produces more z/q, in the L^^^ domain and only reinforces the domain 
structure. 

Therefore, once the adiabaticity condition Eq. (|^) is met, this new channel of sterile 
neutrino production may be potent enough to bring the sterile neutrinos into equilibrium 
with active neutrinos. (In fact, this condition is conservative because neutrinos may cross 
multiple domain boundaries within a Hubble time.) On the other hand, the observationally- 
inferred primordial ^He abundance [lllJl^ and deuterium abundance |13| constrain the total 
number of neutrino flavors in equilibium to be < 3.3 [p|,|6|,p!^ {N^, represents relativistic 
degrees of freedom in neutral fermions) . This constraint implies that the new sterile neutrino 
production channel cannot be efficient before the decoupling of the Ua^a pair production 
process. For a = fi and r, this decoupling temperature is ~ 5 MeV. A constraint on the 
two-family z/^.^ ^ mixing can therefore be obtained by requiring that the adiabaticity 
condition Eq. is not satisfied at T ~ max(5, |5m^/4eV^|^/^) MeV (the latter term in the 
bracket is the temperature at which the growth of L^^ stops): 

\6m'^\ sin^ 29 < 7 x IQ-^eV^ for \6m'^\ < 2.5 x lO^eV^; 

II I l~ 

sm'^2e < 3 X for \6m'^\ > 2.5 x lO^eVl 

For a = e, the constraint from BBN is more severe. The sterile neutrino production 
cannot be efficient not only above the z/gZ/e pair production decoupling temperature T ~ 3 
MeV, but also at the weak freeze-out temperature of T ~ 1 MeV. At this temperature, a 
significant —>■ Vg and Vf. — > Vg transition would cause a deficit in the z/gZ/e number density, 
which cannot be replenished by pair production. A significant deficit in the z^e^'e number 
density causes the neutron-to-proton ratio to freeze out too early and results in a primordial 
^He abundance that is too large. (For example, a 10% deficit in the v^Vf, number density has 
roughly the same effect on the ^He abundance as 3.5.) Therefore the BBN constraint 

on the two-family z/g Vg mixing is obtained at T ~ max(l, |5m^/4eV^|^''^) MeV: 
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\5m'^\sm^2e < 5 x IQ-^eV^ for \5m'^\ <4eV2; 

(10) 

&iYi^2d < 10-^ for \5ni^\ > 4eVl 

These bounds are summarized in Figure 1. They apply in addition to the previous bounds 
based on a universe with homogeneous lepton numbers, and together they offer much tighter 
constraints on the two-family active-sterile neutrino mixing. 

Intriguing results can also be obtained if there is active-sterile neutrino mixing involving 
three or more families. One example is the proposal that a, Vs and z/^ ^ Vs' (in 

principle Vgi and Vs can be the same flavor) mixing might be able to simultaneously explain 
the Super Kamiokande atmospheric neutrino data and satisfy the BBN bound @J^. A stand- 
alone ^ Ug oscillation solution to the Super Kamiokande data would violate the BBN 
bound by bringing z/^ into equilibrium during BBN The double active-sterile neutrino 
oscillation proposal argues that this violation of BBN bound may be avoided if a resonant 
Ut- <-> Ug' transformation in the early universe generates a L^^ that hinders the Ug production 
from the ^ Vg mixing (by creating a Ly^ term in Eq. p|). This argument is no longer 
valid once we consider the existence of the L^^ domains as a result of the resonant Vgi 
transformation. Rather the contrary is true. The Ly^ domains facilitate the production of 
Vg via resonant Vg transformation at domain boundaries. The adiabaticity condition, 

Eq. (^, is modified in this double mixing situation to be: 



2 , „ ./EV . / T 
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where 5m\ ^ |m^^ — 9i is the ^ Vg vacuum mixing angle, and 5m2 ~ m^^ — m^ ^. 
To be consistent with BBN, the adiabaticity condition cannot be satisfied for the double 
mixing system from the onset of the L^^ generation T^. ^ 22\5m\/leV'^\^^^ MeV to the v^^t 
decoupling temperature T ~ 5 MeV. However, for 5m\ ~ 10^^ to 10~^ eV^ and sin^ 26i ~ 1 
(the parameters required to explain the Super Kamiokande data), the adiabaticity condition 
is always satisfied in this double mixing proposal for any reasonable choices of the tau 
neutrino mass. Therefore, BBN unambiguously rules out an active-sterile neutrino oscillation 
explanation to the Super Kamiokande data. 
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Another interesting situation involving multi-family active-sterile neutrino mixing arises 
from neutrino mixing schemes proposed to explain simultaneously the Los Alamos Liquid 
Scintillator Neutrino Detector (LSND) Ve signal, Super Kamiokande atmospheric deficit, 
and solar neutrino deficit ||T5|. In these models, ^ Ve mixing (with m^^ ~ ^-^ ^'^ 

eV^ and sin^ ~ 10~^) is employed to explain the LSND result and ^ mixing (with 
f^l - fnl ~ 10"^ eV^ and sin^26'eo ~ 10"^ for the MSW solution, and \ml - mfj ~ 10"^° 
eV^ and sin^ 20es ~ 1 for the vacuum solution) is invoked to explain the solar neutrino data. 
If there is mixing between z/^ and Us as well, however, with sin^ 2^^^ ^ 10^^^ HjH) the L^^ 
background and domains generated by the mixing in the early universe would imply that 
the Vf. ^ Ug mixing would not populate enough Ug to violate BBN constraints only if 



2 2 

m,, — m,, 



\m^29es < 2 X 10"^^ 



2 2 



1/4 

eV^ ~ 2 X IQ-^^eVl (12) 



leV^ 

(This is in analogy to the previous example if we take 5m\ ^ |m^^ — m'^J and Sm"^ 
\Tnl^ — mlj in Eq. [|lT|.) This requirement is not satisfied by the MSW u^. ^ Vg solution to 
the solar neutrino problem. (Note that if the solar neutrino problem is solved by a vacuum 
Uf. ^ Ug mixing, with |m^^ — m^J ~ lO^"*^" eV^, Eq. (^) will be satisfied and the z/g Ug 
resonant transition will not be adiabatic.) 

Therefore, in light of LSND, Super Kamiokande and solar neutrino experiments, the 
neutrino oscillation explanations of the LSND data and the MSW solution to solar neutrino 
data is inconsistent with BBN unless the z/^ Vg mixing is extremely small, sin^ 2^^^ < 
10^^^. This results holds despite the possibility that the ^ Ve mixing amplitude and the 

^ mixing amplitude could be > 10* times larger. This restriction severly constrains 
the VT--v^-Ve-i^g mixing matrix required to fit the current neutrino experiment results. 

In summary, we have discussed the existence of leptonic domains as an inevitable con- 
sequence of resonant active-sterile neutrino oscillation mechanisms for generation of lepton 
number. Resonant MSW conversion due to the lepton number gradients at domain bound- 
aries therefore provides a new channel for sterile neutrino production. As a result, the 
Big Bang Nucleosynthesis constraint on active-sterile neutrino mixing becomes much more 
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stringent. Likewise for the constraint on multi-family neutrino mixing schemes involving 
sterile neutrinos. We have found that the Ug explanation to the Super Kamiokande 

data is inconsistent with Big Bang Nucleosynthesis in spite of lepton number asymmetries 
generated by other active-sterile neutrino oscillations. We have also found that together 
the Ufj^ <-> explanation of the LSND result and the MSW ^ Ug solution to the solar 
neutrino problem are incompatible with Big Bang Nucleosynthesis considerations unless the 
amplitude of the mixing between i/^ and Us is > 10^ smaller than that between i/f^ and Ue 
and that between and Vs- 
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Figure Captions: 

Figure 1. Parameter spaces to the right of the hatched lines are excluded by BBN. The 
solid lines indicate bounds obtained in this work, and the dashed lines are previous bounds 
assuming a universe with a homogeneous lepton number 
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